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ABSTRACT 



Studies were made of ASCA spectra of seven ultra-luminous compact X-ray 
sources (ULXs) in nearby spiral galaxies; M33 X-8 (Takano et al. 1994), M81 
X-6 (Fabbiano 1988b; Kohmura et al. 1994; Uno 1997), IC 342 Source 1 (Okada 
et al. 1998), Dwingeloo 1 X-l (Reynolds et al. 1997), NGC 1313 Source B 
(Fabbiano & Trinchieri 1987; Petre et al. 1994), and two sources in NGC 4565 
(Mizuno et al. 1999). With the 0.5-10 keV luminosities in the range 10 39-40 
ergs s -1 , they are thought to represent a class of enigmatic X-ray sources 
often found in spiral galaxies. For some of them, the ASCA data are newly 
processed, or the published spectra are reanalyzed. For others, the published 
results are quoted. The ASCA spectra of all these seven sources have been 
described successfully with so called multi-color disk blackbody (MCD) emission 
arising from optically-thick standard accretion disks around black holes. Except 
the case of M33 X-8, the spectra do not exhibit hard tails. For the source 
luminosities not to exceed the Eddington limits, the black holes are inferred 
to have rather high masses, up to ~ 100 solar masses. However, the observed 
innermost disk temperatures of these objects, T in = 1.1 — 1.8 keV, are too 
high to be compatible with the required high black-hole masses, as long as 
the standard accretion disks around Schwarzschild black holes are assumed. 
Similarly high disk temperatures are also observed from two Galactic transients 
with superluminal motions, GRO 1655-40 and GRS 1915+105. The issue of 
unusually high disk temperature may be explained by the black hole rotation, 
which makes the disk get closer to the black hole, and hence hotter. 



Subject headings: black hole physics — galaxies: spiral — X-rays: galaxies 
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1. Introduction 

X-ray radiation from a spiral galaxy consists of emission from discrete X-ray sources, 
hot gaseous media, and in some cases the active galactic nucleus (AGN); for a review, 
see Fabbiano (1988b). Majority of the discrete X-ray sources are presumably accreting 
collapsed objects, in particular low-mass X-ray binaries (LMXBs; a close binary consisting 
of a weakly- magnetized neutron star and a low-mass star). In fact, the X-ray emission 
from M31, with the total 2-10 keV luminosity of 5 x 10 39 ergs s _1 , is dominated by some 
100 LMXBs (Van Speybroeck et al. 1979; Fabbiano, Trinchieri & van Speybroeck 1987; 
Makishima et al. 1989). The most luminous ones in M31 have X-ray luminosities close 
to the Eddington limit for a 1.4 M & neutron star, L^ s = 2 x 10 38 ergs s _1 (Supper et al. 
1997). Similarly, discrete X-ray sources in M33 have X-ray luminosities below L^ s (Long et 
al. 1996), except the one (called X-8) to be discussed later. 

In many other nearby spiral galaxies, however, there reside point-like off-center X-ray 
sources of which the X-ray luminosity significantly exceeds L^ s (e.g., Long 1982; Fabbiano 
1988b; Marston et al. 1995; Colbert & Mushotzky 1999). Actually, the X-ray luminosity 
function of point-like non-AGN sources in nearby spiral galaxies derived with ROSAT 
extends well beyond L E and reaches even ~ 2 x 10 40 ergs s _1 (Read et al. 1997). Although 
some of these luminous X-ray objects are identified with young supernova remnants (e.g., 
Fabbiano & Trinchieri 1987; Schlegel 1994), others have not been identified securely in 
other wavelengths. These luminous non-AGN and non-supernova-remnant X-ray objects 
may be called "ultra- luminous compact X-ray sources" (hereafter abbreviated as ULXs). 
The ULXs cannot generally be explained as collections of multiple sources with luminosity 
below L^ s each, since many of the ULXs exhibit significant time variability. Statistical 
calculations show that these objects are not likely to be less luminous foreground objects or 
background AGNs, either (Fabbiano 1988b). 

Assuming that ULXs are single compact objects powered by mass accretion, there 
are basically two alternative ways of explaining them. One is to regard a ULX as a 
sub-Eddington binary involving a relatively massive black hole (BH) with mass up to 
~ 100 M Q , because the Eddington limit for a body of mass M is given, when electron 
scattering dominates the opacity, as 

L E = 1.5 x 1O 38 (M/M ) ergs s" 1 . (1) 

Here we assume a spherical symmetry, and the hydrogen to helium composition ratio of 
0.76:0.23 by weight. The BH interpretation may be supported by a signature of ULXs as a 
young population; they are more abundant in spirals with later morphological types, and 
they often locate in spiral arms or HII regions (Fabbiano 1988b). However, there is not 
necessarily a consensus on the presence of 100 M Q BHs. 

The other is to regard a ULX as a binary where the X-ray emission is highly collimated 
toward us, or a super-Eddington radiation is maintained. Then, the object can be accreting 
BHs in an ordinary mass range, or even LMXBs. This interpretation has some empirical 
support, because LMC X-2, an LMXB in the Large Magellanic Cloud, exhibits a marginally 
super-Eddington luminosity of 2 x 10 38 ergs s -1 (Tanaka 1997). However, the difficulty of 
this alternative is that there is no known mechanism of producing such an efficient X-ray 
beaming, or a steady super-Eddington radiation. 

Thus, the nature of ULXs has remained a big mystery in the modern X-ray astrophysics. 
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It affects our basic understanding of X-ray properties of normal spiral galaxies, because a 
single ULX could overwhelm the summed contribution from the entire ordinary LMXBs in 
the galaxy. Furthermore, a ULX located near the galaxy center would (and actually does; 
Mizuno et al. 1999; Colbert & Mushotzky 1999) mimic a low-luminosity AGN. 

Clearly, a key to the issue is provided by detailed X-ray spectroscopy, which has 
been enabled by ASCA (Tanaka, Inoue, & Holt 1994). In § 2, we briefly summarize the 
"multi-color disk" (MCD) modeling of emission from optically-thick accretion disks around 
BHs. We show in § 3 that this particular spectral model gives adequate description of 
the ASCA spectra of a fair number of ULXs, on condition that the disk temperature is 
made rather high up to ~ 2 keV. We show in § 4 that the high disk temperature is also 
observed from Galactic BHB binaries (BHBs) with superluminal jets. In § 5, we point out 
a particular problem associated with the BH interpretation of ULXs; that they exhibit 
disk temperatures which are too high for the large BH masses required by their high 
X-ray luminosities. We present various attempts to solve the problem, finally arriving at a 
hypothesis of spinning stellar-mass BHs. Throughout the paper, errors and uncertainties 
refer to 90% confidence limits unless otherwise stated. 



2. The Multicolor-Disk (MCD) Formalism of Accretion Disk Emission 

2.1. Assumed spectral model 

The X-ray emission from a BHB, in so called soft (or high) state, consists of a bright 
soft component (White & Marshall 1984) with a characteristic temperature of ~ 1 keV, and 
a hard tail. The soft component is thought to originate from an optically thick accretion 
disk around the BH, of which the standard accretion disk model by Shakura & Sunyaev 
(1973) provides a prototypical formalism. There are in fact many recent progresses in the 
theory of accretion disks, including the concept of slim accretion disk (Abramowicz et al. 
1988) and advection-dominated accretion flow (ADAF; Narayan & Yi 1995). A unified 
picture is presented by by Esin, McClintock, & Narayan (1997). Nevertheless, as long as 
the high (soft) state of BHBs are concerned, there is no strong observational evidences 
for deviations from the standard Shakura & Sunyaev type accretion disk. Therefore, we 
take the standard accretion disk as a start point, and employ its simple mathematical 
approximation called multi-color disk blackbody (MCD) model (Mitsuda et al. 1984). The 
model is known to give a successful physical description of the X-ray spectra of many BHBs 
(Makishima et al. 1986; Ebisawa et al. 1993; Mineshige et al. 1994; Tanaka & Lewin 1995; 
Tanaka & Shibazaki 1996; Zhang, Cui & Chen 1997; Dotani et al. 1997; Kubota et al. 
1998). 

The MCD model is a superposition of many blackbody elements, up to a certain 
maximum color temperature (sometimes simply called dsik temperature) Tj n that is 
expected to occur near the innermost disk boundary. The local disk temperature T(R) is 
considered to scale as 

T(R)ocR^ 4 , (2) 

where R is radial distance from the BH. Although an exact treatment of the disk emission 
must take into account the inner disk boundary condition and relativistic effects (e.g. 
Ebisawa, Mitsuda & Hanawa 1991), the MCD formalism is known to give a reasonable 



- 5 - 



approximation (e.g., Ebisawa 1991; Dotani et al. 1997; Kubota et al. 1998), providing 
a simple physical insight and a straightforward comparison among different BHs. In 
particular, the validity of equation (fj) has been confirmed observationally in the case of 
Nova Muscae (Mineshige et al. 1994). 



2.2. Assumed geometry 

Assuming a flat disk geometry with inclination i and distance D, the bolometric 
luminosity of an optically thick accretion disk is given as 

L bol = 2-kD 2 ■ f hol (cosi) -1 , (3) 

where /boi is the bolometric flux calculated from the observed band-limited flux, via 
bolometric correction using the MCD model. This L^ Q \ is related to the maximum disk 
color temperature T in and the innermost disk radius _R in as 



^bol 



47r(i? in /OM7W*0 4 • (4) 



Here a is the Stefan-Boltzmann constant, k ~ 1.7 (e.g. Shimura & Takahara 1995) is 
ratio of the color temperature to the effective temperature, or "spectral hardening factor" , 
and £ is a correction factor reflecting the fact that T in occurs at a radius somewhat larger 
than R m . Kubota et al. (1998) give f = (3/7) 1 / 2 (6/7) 3 = 0.412. Because T in is directly 
observable and L bol can be estimated through equation ([|), we can solve equation ([|) as 



D t 2 / -^bol t 2 D / /bol / K s 

Rin = J ^4 = ^ ■ -== ■ J — -j- , 5 

which provides us with the basic tool to estimate the innermost disk radius R m . 

When R in is thus estimated, we may identify it with the radius of the last stable 
Keplerian orbit (e.g. Makishima et al. 1986; Tanaka & Shibazaki 1996). For a non-spinning 
BH of mass M, this condition takes place at 3 R$, where 

2 GM ( M \ , 

fls = _ = 2.95(-) km (6) 

is the Schwarzschild radius, G is the gravitational constant, and c is the light speed. Thus, 
we may in general write as 

.Rin = 3aR s = 8.86a [^~ \ km ( 7 ) 

using a positive parameter a, with a = 1 corresponding to the Schwarzschild BH. 
Conversely, we may obtain an "X-ray estimated" BH mass as 

M = M X r = Rm /8.86a M . (8) 
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2.3. Temperature-luminosity relations 

By substituting equation ([?]) into equation (|), we obtain a temperature-luminosity 
relation for an accretion disk around a BH of mass M, as 



'bol 



^38 / £ 1 ( K \ 4 „2 ( M \ f Ti n \ i 



7.2 x 10 J8 -M — « — — f^) erg s" 1 . (9) 
V 0.41 / VI. 7/ V 10 VkeVy & v; 



Note that the mass accretion rate is unspecified in this relation. 
It may be convenient here to write L\ M \ as 

L hol = r]L E oc V M , (10) 

where a non-dimensional parameter rj means the disk bolometric luminosity normalized 
to the Eddington luminosity of equation (jl]). Using this and equation ([I]), we can then 
eliminate M from equation (^), to obtain another form of temperature-luminosity relation 

as 

^oi = 3.1xl0 39 f4rVf^;V«-Vf7%V 4 ergs- 1 . (11) 



v 0.4iy VI. 7 J ' VkeV/ 
Note that the BH mass is unspecified in this relation. 

By equating equations (|9]) and (|TTD, and solving it for T in , we obtain an important 
scaling of 

Tta = 1 . 2 ( JJ\"(JL) a -^v<(«Y l,t keV . (12) 

Thus, a heavier BH tends to show a lower disk temperature. 



3. Analysis of the ULX Spectra 

Although ASCA has a rather limited angular resolution, we can utilize its wide-band 
(0.5-10 keV) spectroscopic capability for those ULXs which are relatively isolated from 
other X-ray sources in the target galaxies, particularly when helped with the ROSAT 
and/or Einstein images. The ULXs so far studied with ASCA include; the central source 
(X-8) in M33 (Takano et al. 1994), so called Source 1 in IC 342 (Okada et al. 1998), source 
X-6 in M81 (Kohmura 1994; Uno 1997), two sources (A and B) in NGC 1313 (Petre et al. 
1994), one (called X-l) in the newly discovered spiral galaxy Dwingeloo 1 (Reynolds et al. 
1997), and two ULXs in NGC 4565 (Mizuno et al. 1999). Below we review these results, 
and in some cases, reanalyze the published spectra or newly process the ASCA data. In 
Table [l], we give a brief summary of these results. 



EDITOR: PLACE TABLE | HERE. 
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3.1. M33 X-8 

The source X-8 (Long et al. 1981; Gottwald, Pietsch & Hasinger 1987; Trinchieri, 
Fabbiano & Peres 1988) at the center of M33 has a 2-10 keV luminosity reaching 1 x 10 39 
ergs s -1 , which makes it the most luminous X-ray source in the Local Group. In spite of the 
positional coincidence with the optical nucleus, X-8 is unlikely to be a low-lumiosity AGN, 
because of the 106 day periodicity found in its i2(3SViT lightcurve (Dubus et al. 1997), the 
lack of nuclear activity in M33 in other wavelengths, and the stellar kinematic evidence for 
the absence of a massive BH (Lauer et al. 1998). Then, considering its luminosity exceeding 
L^ s , X-8 may be interpreted alternatively as a ULX. 

Takano et al. (1994), and subsequently Mizuno (2000), analyzed the spectra of M33 
X-8, taken with the Gas Imaging Spectrometer (GIS; Ohashi et al. 1996; Makishima et 
al. 1996) and the Solid-State Imaging Spectrometer (SIS; Burke et al. 1994; Yamashita 
et al. 1997) onboard ASCA. The spectra have been described successfully with an MCD 
model of T in ~ 1.15 keV, as quoted in Table |l], on condition that an additional power-law 
component is included to describe a hard tail above ~ 5 keV. Although the MCD emission 
is associated with either LMXBs or BHBs, X-8 is too luminous for an ordinary LMXB, and 
its power-law hard tail is characteristic of BHBs. Takano et al. (1994) and Mizuno (2000) 
therefore conclude that X-8 is a mass-accreting stellar-mass (~ 10 M ) BH. These results 
have been reconfirmed by Colbert & Mushotzky (1999). 

Quantitatively, if we adopt i = 0°, the measured values of T in = 1.15 keV and 
Lboi = 3.8 x 10 38 ergs s -1 of X-8 (Table [I]) can be consistently reproduced by substituting 
M = 5.5 M & and t] = 0.46 into equations ([|) and (|ll|), together with canonical values 
of a = 1, £ = 0.41, and k = 1.7. If instead, e.g., i = 60° is adopted, L bol is doubled to 
7.6 x 10 38 ergs s _1 , and the solution becomes M ~ 7.8 M & and i] ~ 0.65. These estimates 
reconfirm the inference made by Takano et al. (1994). We are then inspired to regard M33 
X-8 as the nearest and prototypical ULX. 



3.2. IC 342 Source 1 

Another case of significant interest is "Source 1" (Fabbiano & Trinchieri 1987; Bregman 
et al. 1993; Makishima 1994) located in a spiral arm of the galaxy IC 342. As reported by 
Okada et al. (1998), its spectrum obtained with ASCA can be expressed by an MCD model 
of maximum color temperature T- m ~ 1.8 keV, whereas the power-law and Bremsstrahlung 
models disagree with the data. We here reanalyze the same ASCA data in the same way 
as Okada et al. (1998). Detailed study of the time variation is reported elsewhere (Mizuno 
2000). 

The derived SIS and GIS spectra of Source 1 are presented in Figures [p. Here and 
hereafter, all the ASCA spectra are presented after background subtraction, but without 
removing the instrumental responses. We fitted the SIS/GIS spectra jointly, using three 
different spectral models; power-law, thermal Bremsstrahlung, or MCD, all modified at low 
energies with photoelectric absorption. Then, as shown in Table |2|, only the MCD model 
among the three has given an acceptable joint fit, in agreement with Okada et al. (1998). 
Predictions of the best-fit power-law and MCD models are shown respectively in Figure [l]a 
and Figure 0b as solid histograms. 
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The successful MCD fit encourages us to interpret this source as an accreting BH, like 
M33 X-8. However, the BH must be fairly massive (Table [3]) for the high source luminosity 
to stay within the Eddington limit. More importantly, the measured disk temperature, 
T in ~ 1-8 keV, much exceeds those (0.4-1.0 keV) typically found among established BHBs 
in the Milky Way and Magellanic clouds, in contradiction to the prediction of equation fli~2"D 
that more massive BHs should have lower disk temperatures. These problems, pointed out 
by Okada et al. (1998), Mizuno et al. (1999) and Mizuno (2000), are found commonly with 
other ULXs, as revealed hereafter. 

We further tried two other spectral models of convex shapes. One is unsaturated 
Comptonization (UC) model by Sunyaev & Titarchuk (1980), describing thermal electrons 
of temperature T e Compton up-scattering soft seed photons into X-rays. It has given 
a fairly good fit to the spectra of IC 342 Source 1, although not as good as the MCD 
model (Table [J). This is of no surprise, because the MCD and UC models can take nearly 
identical shapes (Ebisawa 1991), when T e is made close to T in and the optical depth for 
scattering, r es , is adjusted. The other is a broken power- law, often used to approximate 
synchrotron radiation spectra emitted from jet-dominated sources. This model, modified 
with photoelectric absorption, gave a value of chi-square in between those from the MCD 
and Bremsstrahlung fits (Table |2|). Physical meanings of these alternative modelings are 
examined in § 5.1. 

EDITOR: PLACE FIGURE |T] HERE. 
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3.3. M81 X-6 

Because of the explosion of the supernova SN1993J, the region around the spiral 
galaxy M81 has been observed repeatedly with ASCA. These observations have provided 
valuable information on SN1993J (Kohmura et al. 1994; Kohmura 1994; Uno 1997), the 
low-luminosity AGN in M81 (Ishisaki et al. 1996; Iyomoto 1999), and the previously known 
ULX called M81 X-6 (Fabbiano 1988a) which lies close to SN1993J. A merit of X-6 is the 
accurate knowledge of the distance to the host galaxy, 3.6 ± 0.3 Mpc (Freedman et al. 
1994). 

Figure [T|c presents the SIS spectrum of X-6 (plus SN1993J), produced by Uno (1997) 
who screened the data through standard procedure, accumulated photons falling within 1/5 
of X-6, and then subtracted background as well as contamination from the M81 nucleus. 
The spectrum co-adds two observations, conducted on 1994 April 1 through 2 for 28 ks 
and on 1994 October 21 through 22 for 37 ks. The supernova had faded away significantly 
by these observations, and its contribution to the present spectrum can be expressed by a 
steep power law having a photon index of 3.0 and a photon flux density of 1.0 x 10~ 4 ergs 
s _1 cm -2 keV -1 at 1 keV (Uno 1997). Therefore, the contamination is limited to energies 
below ~ 2 keV. In the present study of M81 X-6 we utilize this SIS spectrum only, since the 
poorer spatial resolution of the GIS increases the supernova contamination. Further data 
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analysis of X-6, including the intensity-sorted spectroscopy, is reported elsewhere (Mizuno 
2000). 

We fitted the SIS spectrum with five alternative spectral models; power-law, thermal 
Bremsstrahlung, MCD, UC, or broken power-law, all modified at low energies with 
photoelectric absorption. The SN1993J contribution was taken into account as a fixed 
model component mentioned above. As summarized in Table ^, the results have been 
qualitatively similar to those from IC 342 Source 1: the power-law model failed, the 
Bremsstrahlung fit is much better but still unacceptable, while the remaining three models 
have been generally successful. In any case, reservation should be put on the obtained N^, 
because of the supernova contamination below ~ 2 keV. 

Based on the MCD fit and equation the disk bolometric luminosity becomes as 
given in Table [1]; clearly, X-8 is a ULX, and thanks to the accurate distance, its luminosity 
is considerably more reliable than that of IC 342 Source 1. The BH mass required by the 
Eddington limit (Table [I]) may be reasonable as long as i is rather small. However, in 
reference to equation (|l^), the measured value of T- m (Table is too high for the rather 
high BH mass. 

The high value of T in suggests the presence of a separate harder spectral component, 
e.g., a power-law hard tail normally seen from soft-state BHBs just as is the case with M33 
X-8 (§ 3.1). Accordingly, we refitted the X-6 spectrum by a sum of an MCD continuum, 
and a hard power-law of which the photon index is fixed at 2.3, typical of high-state BHBs. 
However the data did not require the hard component, with its 0.5-10 keV flux being at 
most 18% of the total source flux in the same range. Furthermore, inclusion of the hard 
tail component at the allowed upper limit increased the MCD temperature slightly, from 
1.48 keV to 1.52 keV. Alternatively, a hard blackbody (BB) component of temperature ~ 2 
keV might be present, like in X-ray luminous LMXBs (Mitsuda et al. 1984; Tanaka 1997). 
Accordingly, we replaced the power-law hard tail with a blackbody of temperature fixed at 
2.0 keV. However, the data did not require the BB component either, and its contribution 
to the overall 0.5-10 keV flux turned out to be at most 14%. 



3.4. Dwingeloo 1 X-l 

Through ASCA observations of the recently discovered nearby spiral galaxy Dwingeloo 1 
located behind the Milky Way, Reynolds et al. (1997) detected a ULX, and named it X-l. 
Its spectrum, presented in Figure |I]d, was described well with several alternative models 
because of rather poor data statistics, but the MCD fit was not attempted by Reynolds et 
al. (1997). We hence reanalyze the same GIS and SIS spectra, in terms of the MCD model 
as well as power-law and Bremsstrahlung models. As shown in Table |2|, all these models 
have been acceptable, with the power-law and Bremsstrahlung results being consistent with 
those reported by Reynolds et al. (1997). The MCD model is somewhat preferred among 
the three, because the absorption associated with it agrees with the line-of-sight Galactic 
value of A% = 7 x 10 21 cm -1 (see Reynolds et al. 1997); the other two models require 
excess absorption which could be artificial. The source again exhibits a quite high disk 
temperature, ~ 1.8 keV. 
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3.5. NGC 1313 Source B 

The spiral galaxy NGC 1313 is known to contain three bright X-ray objects; one 
identified with SN1978k, another called Source A located near the nucleus, and the other 
called Source B associated with a spiral arm (Fabbiano & Trinchieri 1987). Here we focus 
on Source B which is regarded as a ULX (Petre et al. 1994; Colbert et al. 1995); Source A 
could be composite. The NGC 1313 region was observed twice with ASCA, on 1993 July 
12 and 1995 November 29. 

The first ASCA observation was reported by Petre et al. (1994), who found that 
the Source B spectra can be described well with a Bremsstrahlung or a Raymond-Smith 
model (Raymond & Smith 1977), whereas the power-law model is rejected. However, 
they did not try an MCD fit. We therefore reanalyzed the same SIS and GIS spectra of 
NGC 1313 Source B, which are reproduced in our Figure [l|e. As summarized in Table 0, 
the power-law fit is unacceptable, and the Bremsstrahlung fit is marginally acceptable; the 
derived parameters agree with those of Petre et al. (1994). The MCD fit is fully acceptable, 
and the obtained disk temperature is again rather high. 

We have also analyzed the ASCA data of Source B from the second ASCA observation, 
by screening the data in the standard way, and accumulating the SIS and GIS photons 
separately over a common region of radius 3' centered on the source. We subtracted the 
GIS background using blank-sky data, and the SIS background using source-free regions of 
the same on-source data. As shown in Figure [l|f and Table ^, the MCD fit is again most 
successful. As listed in Table |2], the source flux decreased by a factor of 2.5 from 1993 to 
1995, meanwhile the disk temperature decreased by a factor of ~ 1.4. There is marginal 
evidence (Table |l]) for an increase in the disk radius as the source became fainter. This 
phenomenon is reported in further detail by Mizuno (2000). 



3.6. Two ULXs in NGC 4565 

Finally, we refer to the recent results by Mizuno et al. (1999), who studied X-ray 
emission from the edge-on spiral galaxy NGC 4565 using ASCA and ROSAT. The galaxy 
hosts two bright point-like X-ray sources, the fainter one coincident with the nucleus 
(hereafter called center source) while the brighter one displaced ~ 2 kpc above the galaxy 
disk (hereafter off-center source), both showing high bolometric luminosities (Table 0). 
Their ASCA spectra have been described successfully by the MCD model with rather high 
disk temperatures (Table 0), and marginally well by the Bremsstrahlung model, but the 
power-law fit was unacceptable for the off-center source. Based on the absence of significant 
spectral absorption, Mizuno et al. (1999) conclude that the center source is in fact a ULX 
in NGC 4565 rather than an AGN. 



4. Galactic Transients with Super-luminal Jets 

Although no ULXs are known in our Galaxy or M31, the issue of rather high disk 
temperature is also seen among some Galactic transient sources (Zhang, Cui & Chen 1997). 
The best examples are GRO J1655— 40 and GRS 1915+105, both exhibiting super-luminal 
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radio jets which have allowed accurate determinations of the system inclination and 
distances. We here analyze the ASCA data of these two sources, because their resemblance 
to ULXs is considered to be meaningful. As they are too bright for the SIS, we utilize the 
GIS data only. 



4.1. GRO J1655-40 

The transient source GRO J1655-40 has % = 69.°50 ± 0.°08 (Orosz et al. 1997) and 
D = 3.2 kpc (Hjellming & Rupen 1995). It has a secure BH secondary, with the optically 
determined BH mass M opt of 7.0 ± 0.2 M (Orosz et al. 1997) or 5.5 - 7.9 M (Shahbaz et 
al. 1999). The RXTE observations are reported by Mendez et al. (1998). 

We reanalyzed the deadtime-corrected ASCA GIS spectrum of this source, acquired on 
1995 August 15-16 and published by Ueda et al. (1998) as their Figure 2a; they fitted it 
by a power-law model with a high-energy cutoff, to focus upon the iron absorption feature. 
We fit it with the MCD model, plus a power-law of which the photon index is fixed at 2.36 
after Zhang et al. (1997). The iron absorption feature was represented with a spectral 
notch at 6.95 keV according to Ueda et al. (1998). We have obtained an acceptable fit with 
X 2 '/v = 91.6/91, in spite of very high signal statistics. This reconfirms the correctness of 
the spectral modeling, and the accuracy of the instrumental calibration (Makishima et al. 
1996). The derived MCD parameters are given in Table [lj the observed disk bolometric 
luminosity is ~ 1/8 of the Eddington limit for M opt . The absorption turned out to be 
(9 ± 1) x 10 21 cm" 2 . 

Through equation fl5|), we obtain R m ^/cosi = 11.4 km, which in turn yields an X-ray 
mass estimate of Mxr ~ 2.2 M via equation (|S[). This is in fact much lower than the 
M opt quoted above. This Mxr does not change by more than 10% even when forcing the 
power-law normalization to be zero (then the fit becomes unacceptable), or leaving the 
power-law slope free. Thus, the disk temperature is apparently too high, and the disk 
radius is too small, just like in the case of the ULXs described in § 3. 

We also attempted the spectral fitting using general relativistic accretion disk (GRAD) 
model (Hanawa 1989; Ebisawa, Mitsuda Sz Hanawa 1991), which improves over the MCD 
formalism by taking into account the inner boundary condition and relativistic effects. It 
has the BH mass M and the accretion rate M as free parameters, with i and k adjustable 
auxiliary parameters. We fixed i at 69°. 5 and k at 1.7, and again fixed the power-law 
photon index at 2.36. As shown in Figure 0a, an acceptable {x 2 /v — 70.4/91) fit has been 
obtained, with M = 2.91 ± 0.10 M and M = (2.64t£$) x 10 18 g s -1 . The value of N H was 
not much different from that of the MCD fit (Table [I]). Thus, the mass estimate somewhat 
increases, but it still remains significantly below M opt . 

These results reconfirm the report by Zhang, Cui & Chen (1997), who used the ASCA 
data acquired on the same occasion to derive R m = 22 km, employing also k — 1.7 and 
relativistic corrections to the MCD results. This has led Zhang, Cui & Chen (1997) to 
propose that the BH in GRO J1655— 40 is rapidly spinning, and the accretion disk is 
prograde to the BH rotation. We come back to the issue in § 5.6. 
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4.2. GRS 1915+105 

Similarly, the transient GRS 1915+105 is a promising BH candidate, and has been 
studied by a number of authors including Belloni et al. (1997) and Zhang, Cui & Chen 
(1997). It has i = 70° and D = 12.5 kpc (Mirabel & Rodriguez 1994), although the optical 
estimate of the BH mass is not yet available. Belloni et al. (1997) report Lboi ~ 1.4 x 10 39 
ergs s , which indicates Me > 9M . Then, the disk temperature should be < 1.2 keV, 



from equation (|12|) . However, in reality, a much higher value of T in = 2.27 ± 0.03 keV was 



obtained with RXTE during the outburst peak (Belloni et al. 1997). 

We analyzed the ASCA GIS data of GRS 1915+105, acquired on 1995 April 20 for 18 
ks and described partially by Ebisawa (1997) and Kotani et al. (2000). The GIS data were 
screened in the standard way, and the spectra from GIS2 and GIS3 were added together 
after appropriate dead-time corrections (Makishima et al. 1996) as presented in Figure 0b. 
Background is completely negligible. It has been fitted successfully (x 2 '/v = 82.5/81) with 
the MCD continuum, on condition that it is multiplied by the following six factors; (i) a 
photoelectric absorption by = (4.3 ± 0.1) x 10 22 cm -2 ; (ii) a smeared Fe-K edge with 
the center energy 7.12 ± 0.08 keV, optical depth 1.6 ± 0.2, and a fixed width of 10 keV; 
(iii) a narrow K a line of He-like and/or H-like iron, centered at 6.98 ± 0.03 keV with an 
equivalent width (EW) of 55 ± 6 eV, (iv) a narrow Fe Kp line at 8.15 ± 0.09 keV with an 
EW of 20 ± 6 eV, (v) a third narrow absorption line at 1.90 ± 0.02 keV with an EW of 
19 + 2 eV, attributable to K a line of H-like and/or He-like Si; and (vi) a broad absorption 
line at 2.52 ± 0.04 keV, having width of 0.8 keV and an EW of 22 ± 3 eV, presumably a 
blend of Si-Kg and S-K Q lines. No power-law tail was required. 

The derived MCD parameters are given in Table |I|. The bolometric luminosity 
becomes 1.2 x 10 39 ergs s _1 , which is similar to that measured with RXTE. The value of 
Mxr = 3.5 ± 0.4 M (Table [I]) falls significantly below M E . Thus, as reported by Belloni et 
al. (1997), GRS 1915+105 exhibit a very high T in and too small R m , like GRO J1655-40 
and the ULXs. 



EDITOR: PLACE FIGURE | HERE. 



5. Discussion 



5.1. Spectral properties of the sample ULXs 

We have studied the ASCA spectra of seven ULXs, partially referring to published 
results. The Bremsstrahlung model failed for least in two spectra (IC 342 Source 1, and 
M81 X-6), and the power-law fit is generally unacceptable except in the faintest objects. 
In contrast, all the spectra have been described successfully by the MCD model (plus a 
hard power-law in M33 X-8) with a relatively high disk temperature in the range 1.0-1.8 
keV. Thus, the seven ULXs form a rather homogeneous sample, and hence are thought to 
be the same type of objects. With a possible exception of Dwingeloo 1 X-l with rather 
unconstrained spectral modeling, none of them are background AGNs, since no AGN is 
known to exhibit the MCD-type spectra in the ASCA range. 
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We found that the UC model and the broken power-law model can also describe fairly 
well the two ULX spectra of high statistics, those of IC 342 Source 1 (§ 3.2) and M81 X-6 
(§ 3.3). These models will give acceptable fits to the remaining ULX spectra that have 
poorer statistics. However if the UC interpretation were correct, the spectrum should bear 
a prominent Fe-K edge absorption feature at ~ 7 keV, because the large optical depth for 
electron scattering required by the two ULX spectra, r es = 20 — 30, implies a huge hydrogen 
column density of 5 x 10 25 cm -2 , and because the heavy ions would not be completely 
ionized at the suggested electron temperature of T c ~ 1 keV. Evidently, such a feature is 
absent in the observed ULX spectra, making the UC interpretation of the ULX spectra 
physically unrealistic. Similarly, the broken power-law model (appropriate for synchrotron 
emission) is physically unrealistic, since the break in the photon index of 1.1—1.5 required 
by the two sources is too abrupt to be interpreted in terms of synchrotron cooling. 

These considerations strongly suggest that the emission from these ULXs actually 
originates in optically thick accretion disks as approximated by the MCD modeling, and 
hence the objects are either BHBs or LMXBs. The disk temperatures of the ULXs are closer 
to those of luminous LMXBs (Mitsuda et al. 1984; Tanaka 1997), rather than to those of 
BHBs. However, in order for LMXBs to appear as ULXs, the emission must be highly 
anisotropic, and/or super-Eddington by nearly two orders of magnitude. Furthermore, the 
hard BB component characterizing the LMXB spectra is absent, at least in the M81 X-6 
spectrum, with an upper limit of 14% (§ 3.3). This is considerably lower than the fractional 
BB contribution of ~ 50% observed from the most luminous LMXBs (Mitsuda et al. 1984). 
We therefore conclude that the ULXs are accreting BHBs rather than a version of LMXB. 
This makes a very important step towards understanding these enigmatic objects. 



5.2. Problems with the accretion-disk interpretation 

In Table [l], we compile the disk bolometric luminosity Lboi via equation (|3|), the 
minimum BH mass Me necessary to make Lboi stay within the Eddington limit, and the 
disk inner radius i? in given by equation @. The assumed distance D to each host galaxy 
is also given there. For all ULXs, we have assumed the face-on inclination, i.e. i ~ 0, that 
makes Me lowest. This assumption is not necessarily too arbitrary, because face-on systems 
must have higher fluxes for a given bolometric luminosity according to equation (^), and 
hence more selectively detectable, than those with larger inclinations; we hence retain the 
assumption of i ~ for the ULXs throughout the paper. When D and i are changed, these 
quantities scale as 

M E oc /^(cosfT 1 , Mkr oc R m oc D (cosz) -1 ^ 2 . (13) 

In Table [I], Me takes rather high values, reaching 70 — 80 M for IC 342 Source 1 and 
NGC 4565 off-center source. Although there is not yet a consensus on the formation of 
BHs in such a mass range, some (e.g. Fryer 1999) argue that a star heavier than ~ 40 M Q 
can directly form a BH without supernova explosion (and hence without losing much of the 
progenitor mass), and there are a fair number of observations of stars in the mass range of 
100 — 150 M Q (e.g. Krabbe et al. 1995). Therefore, the high BH mass may not be a big 
difficulty 

A more profound problem is that, except M33 X-8 (and NGC 1313 Source B in 1995), 
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the observed T in is too high for the required high BH mass (Okada et al. 1998; Mizuno et 
al. 1999; Mizuno 2000), as mentioned in § 3. The high values are unlikely to be an artifact 
of the data analysis (§ 3), and substitution of these values into equation ([TJ) would make 
the BH masses < 5 M even at the extreme of rj = 1. This makes self-inconsistent the BH 
interpretation of ULXs, where a rather high BH mass is inevitable. This problem is best 
visualized by Figure to be called an X-ray "H-R diagram" for BHs, which summarizes 
various BHBs and ULXs on the plane of T in vs. Lboi- There, a family of grid lines with a 
positive slope (Lboi ^ ^in) represent the loci of constant M, expressed by equation ([]). The 
other family of grid lines with a negative slope (Lboi ^ ^ln 4 ) indicate the loci of constant 
rj, expressed by equation (|TTJ) . In order to consistently interpret an X-ray source as an 
accreting non-spinning BH, the data point must fall on the region of r] < 1 (non- violation of 
the Eddington limit) and M > 3 M (heavier than a neutron star) in Figure [5| The reality, 
however, is that six out of the seven ULXs (except M33 X-8) fall on the super-Eddington 
regime in this plot. As to NGC 1313 Source B, the 1993 data point is super-Eddington, 
whereas the 1995 point is marginally sub-Eddington. 

EDITOR: PLACE FIGURE § HERE. 

The problem may equivalently be stated in the following way. The observed high values 
of T in , together with equation (j^j, yield rather small values of Ri n (Table |I]). Then, the 
BH mass Mxr, determined by equation (H), becomes uncomfortably low. Quantitatively, 
the derived Mxr falls by a factor of 1.6-6.0 short of the minimum BH mass Me to satisfy 
the Eddington limit. The issue may be illustrated in Figure £|, to be called a "mass-radius 
diagram" for BHs, where we compare the values of i?; n against Me or M opt . 

It is important to note that the same problem is found also with the two Galactic jet 
sources studied in § 4. We accordingly plot these two objects as well in Figure |3| (more than 
once for each object) and Figure §. Indeed, in the H-R diagram of Figure [| GRS 1915+105 
falls on the super-Eddington regime, and GRO J1655— 40 is located near the lower-mass 
boundary of a BH. Similarly in the mass-radius diagram of Figure |4|, the two jet sources 
(both referring to our own data) deviate from the relation of ordinary BHBs (see § 5.4), 
and line up better with the ULXs. 

EDITOR: PLACE FIGURE | HERE. 



5.3. Simple solutions 

As first-cut attempts to solve the above problem, let us consider trivial possibilities, 
namely revisions of the inclination i, the Eddington limit, or the distance. 

So far, we have assumed i = 0. However, as is apparent from equation (p~3|) , the ratio 
of Me/Mxr increases with i as oc (cosz)~ 1//2 . Therefore the problem of high Me/Mxr ratio 
gets worse if i is increased from ~ 0. Thus, i has already been selected to make the problem 
least severe. Larger values of i would also make the BH mass still higher. 
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Alternatively, the problems may result from our wrong application of the spherically- 
symmetric Eddington limit of equation (|l|) to the systems considered here, which are 
basically axi-symmetric. However, the critical upper-limit luminosity of an axi-symmetric 
BH system does not differ very much from the spherical one (e.g., Abramowicz et al. 1988; 
Beloborodov 1998). This justifies our use of equation (0). 

Yet another simple possibility is that the source distances are systematically 
overestimated: note Me/Mxr oc D from equation (fl3|). However, in order to bring Me/Mxr 
down to unity, we have to reduce the distances to IC 342, M81, NGC 1313, and NGC 4565 
by a factor of 6.0, 2.2, 3.1, and 3.7, respectively. Such large distance revisions would be 
unwarranted, especially for M81 of which the distance is known to within ±10% (Freedman 
et al. 1994). Therefore, this attempt is unsuccessful, as already noticed by Mizuno et al. 
(1999) regarding the ULXs in NGC 4565. 



5.4. Possible errors in £k 2 

As is obvious from equations (|5|) and (||), the Me vs. Mxr discrepancy may 
alternatively be caused because our choice of £ = 0.41 (Kubota et al. 1998) and k = 1.7 
(Shimura & Takahara 1995), or £k 2 = 1.18, is wrong. In this subsection, we therefore 
attempt to "calibrate" the values of £ and k mainly using established BHBs. 

The factor £ corrects the MCD formalism for the inner boundary condition of the 
accretion disk, in such a way that the disk temperature reaches maximum at -Rm/£ instead 
of R m . We may then examine our choice of £ = 0.41 (Kubota et al. 1998) in reference 
to the GRAD model mentioned in § 4.1 that properly takes into account this effect. We 
accordingly simulated a number of spectra using the GRAD model, by changing M, M 
and i, but fixing k at 1.7. The spectra were analyzed with the MCD formalism, employing 
£ = 0.41 and n = 1.7. Then, over the inclination range of 0° to 75° and over T- m = 1.2 — 2.0 
keV, the X-ray mass Mxr from the MCD analysis agreed within ±25% with the BH mass 
M employed as initial inputs to the GRAD model simulations. 

We can also utilize Cygnus X-l, the prime BHB, for this purpose. Applying the 
GRAD model to the ASCA data of Cygnus X-l acquired in the 1997 May soft state, and 
assuming D = 2.5 kpc and i = 30°, Dotani et al. (1997) derived Mxr = 12±f M . If 
we apply our MCD formalism with C.K 2 = 1.18 to the same dataset for Cygnus X-l, we 
obtain R in = 90 ± 18 km (Table |l|), and hence Mxr = 10 ± 2 M , which is close to the 
GRAD result by Dotani et al. (1997). Furthermore, our value of Mxr is consistent with 
the latest optical estimate of the BH mass in Cygnus X-l, M opt = 10.1±t:| M© (Herreo et 
al. 1995). Thus, our overall knowledge of Cygnus X-l supports the value of £k 2 ~ 1.2. We 
plot Cygnus X-l in Figures || and f|. 

We may further utilize the two Magellanic BHBs, LMC X-l and LMC X-3, which have 
well established distances of 50-55 kpc; here we adopt 55 kpc. We hence analyzed the 
ASCA GIS data of LMC X-l and LMC X-3, obtained on 1995 April 2 for 22 ks and 1995 
April 14 for 23 ks, respectively. Their GIS spectra have both been fitted successfully by a 
model consisting of an MCD component and a power-law tail, whenall model parameters 
are left free to vary. The derived MCD parameters are given in Table III while the power-law 
index turned out to be 2.36±g;i| for LMC X-l and 2.52i£ 2 j for LMC X-3. The photoelectric 
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absorption is JV H = (5.9 ± 0.3) x 10 21 cm- 2 in LMC X-l , and N R < 2.9 x 10 21 cm" 2 
in LMC X-3. The measured disk temperatures, given in Table p], agree with the Ginga 
measurements (Ebisawa 1991; Ebisawa et al. 1993). Using i ~ 60° for LMC X-l (Cowley 
1992; Cowley et al. 1995) and i = 66° ± 2° for LMC X-3 (Kuiper et al. 1988), we obtain 
the disk radii as given in Table [p. In Figure these radii are in good consistency with the 
optical mass estimates of M opt = 4—10 M & for LMC X-l (Cowley 1992; Cowley et al. 
1995) and M opt = 5.0 - 7.2 M for LMC X-3 (Kuiper et al. 1988, D = 55 kpc). 

Thus, the value of £k 2 we have employed is not particularly wrong, and may not be 
changed by more than a factor of 1.5 or so. This tolerance is not large enough to solve the 
Me vs. Mx discrepancy for many of our sources. 



5.5. Luminosity dependences 

Even though we have argued in § 5.4 for the justification of our estimates on R m , the 
disk structure may depend on Lboi or M, particularly when Lboi approaches the Eddington 
limit. The electron scattering effect may get severer, which will increase k, e.g., from 
1.7 to 1.8-2.0 (Shimura & Takahara 1995). Furthermore, the true disk temperature may 
become higher than is described by the standard accretion disk model, because entropy 
of the accreting matter may start being "advected" into the BH as described by the slim 
disk model (Abramowicz et al. 1988). Recently, Watarai et al. (1999) have shown via 
calculation that significant amount of X-rays are radiated even from the region inside 3i? s 
in a slim accretion disk accreting close to the Eddington limit. (We may also refer to 
Reynolds & Begelman (1997) who discuss the effect of accretion flow inside 3R S .) While 
this mechanism can potentially explain what we have observed, its effects will vanish when 
the accretion rate becomes sufficiently low. Therefore, we are urged to search the available 
data for luminosity-dependent effects, such as apparent changes in R [n . 

As mentioned in § 5.4.3, LMC X-3 was observed extensively with Ginga. In Figure 0, 
we therefore plot all the Ginga and ASCA data points for LMC X-3. Thus, the luminosity 
of LMC X-3 varied by a factor of ~ 7, reaching ~ 2/3 of the Eddington limit. However, the 
data points are aligned very well with a grid line for a constant mass as already reported 
previously (Ebisawa 1991; Ebisawa et al. 1993), and the implied mass is consistent with 
M opt . Therefore, we do not see any luminosity-dependent changes in R in . 

The bright BH transient GS 2000+25, discovered with the Ginga ASM (Tsunemi et al. 
1989), is still more suited for the study of the luminosity dependence of accretion disks, 
because of its large luminosity swing. After the discovery, GS 2000+25 was observed with 
the Ginga LAC eight times spanning 240 days, covering almost the outburst peak (Ebisawa 
1991). The obtained 2-30 keV spectra were all expressed well with an MCD component 
plus a hard power-law tail. In Figure |^, we plot seven out of the eight Ginga data points 
(discarding the last one), assuming a source distance of 3 kpc and i = 65° after Callanan 
et al. (1996). Thus, over a very wide range of M up to n ~ 0.7, the data points generally 
distribute along a constant-mass grid line, as already noticed by Ebisawa (1991). The 
value of Mx ~ 7 M indicated by Figure || agrees with M opt measured by various authors 
(Filippenko, Matheson & Barth 1995; Beekman et al. 1996; Callanan et al. 1996), which 
cluster around 5 — 8 M Q . Thus, the standard-disk interpretation is fully self consistent as 
to this objects, with little evidence of luminosity-dependent effects. 
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The superluminal source GRO J1655-40 itself exhibited a large luminosity decline as 
well. In Figure |3|, the ASCA and RXTE data points generally align along a constant-mass 
locus, even though the absolute value of M X r falls considerably below M opt as argued 
repeatedly. In short, the small value of Mxr (or the too high T in ) of this jet source is not 
specific to its high-luminosity states, but appears intrinsic to it. 

We may further argue against luminosity-dependent effects by comparing 
GRO J1655— 40 and LMC X-3. These two objects are very similar in terms of the 
BH mass, inclination, and L^oi at the time of the ASCA observation (Table |l|). Nevertheless, 
their disk temperatures differ systematically by 40%, which makes their disk radii 
significantly different in Figure ||. 

We therefore conclude that luminosity-dependent changes in the accretion disk, if any, 
cannot fully explain the high disk temperature seen among a certain class of BHs. We 
regard it as an intrinsic property of such BHs. 



5.6. The spinning black hole scenario 

We are finally urged to examine the remaining degree of freedom in our MCD scenario, 
i.e., the assumption of a = 1 in equation fl7|), which reflects the fact that the stable circular 
orbit can exist only outside 3 R s in the case of a Schwarzschild (i.e., a non-rotating) BH. 
If instead the BH is rotating significantly (i.e., a Kerr BH), and the accretion disk is 
prograde to it, the last stable Keplerian orbit gets closer to the BH, down to R m = 0.5 R s 
for a maximally rotating BH. As a result, we expect a Kerr BH to have a smaller value of 
a, down to the extreme of a = 1/6, than a Schwarzschild BH of the same mass and the 
same accretion rate. Then, a Kerr BH can attain a higher value of T in , neglecting for the 
moment stronger relativistic effects. In Figure |], the dotted line shows the prediction for 
the maximally rotating BHs, obtained by substituting a = 1/6 into equation (0). Thus, 
except for IC 342 Source 1 for which the distance is rather uncertain, all the objects now 
fall on the physically reasonable region. Therefore we propose that the ULXs contain Kerr 
BHs. This idea was first proposed by Zhang, Cui & Chen (1997), to explain the high disk 
temperatures of GRS 1915+105 and GRO J1655— 40, and subsequently applied by Mizuno 
et al. (1999) to the two ULXs in NGC 4565. 

Of course, investigation of the emission from rotating BHs must fully take into account 
relativistic effects. Such studies have been conduced by several authors, including, e.g., 
Asaoka (1987), Zhang, Cui & Chen (1997), and Beloborodov (1998). Asaoka (1987) showed 
that the spectrum from an optically-thick accretion disk around a Kerr BH is similar in 
shape to that around a Schwarzschild BH, except a higher temperature. Zhang, Cui & 
Chen (1997) calculated how the maximum color temperature of the accretion disk, T- m , 
depends on the BH rotation. They show that, when the disk is prograde to a maximally 
rotating BH, T in can be higher by a factor of ~ 3 than that around a Schwarzschild BH 
with the same BH mass and the same mass accretion rate. This is large enough to account 
for the high disk temperature of our sample objects. (This temperature enhancement factor 
becomes 6 0,75 = 3.8 in terms of the MCD formalism which neglects the general relativistic 
effects.) 

There is another merit in appealing to the BH rotation, that the radiation efficiency, 
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defined as the ratio L^ \/Mc 2 , gets higher as the BH angular momentum increases. The 
efficiency is 0.057 for a Schwarzschild BH, while it increases to ~ 0.4 for a maximally 
rotating BH fed with a prograde disk (see Zhang, Cui & Chen 1997). Therefore, in 
a maximally rotating BH, the same mass accretion rate can produces 7 times higher 
luminosity than in a non-rotating BH, although the same Eddington limit should still apply. 

How to make rather massive (several tens M & ) Kerr BHs is an intriguing issue. There 
seems at least three different scenarios. One is the simplest scenario, to assume that a 
massive star with a rapid rotation collapses into a spinning BH. Another is to invoke a 
merger of two ordinary stellar-mass BHs, where the angular momentum is brought in from 
their orbital motion, and the BH mass is doubled. The other is to assume that a slowly 
spinning BH is spun up due to angular momentum supplied by the accreting matter, just 
as binary X-ray pulsars are spun up. Among the three scenarios, the last one is particularly 
attractive because it can potentially explain the large BH mass and the high angular 
momentum at the same time. Although Chen, Cui & Zhang (1997) argued that this 
mechanism is inefficient, the scenario may change significantly if the accretion is allowed to 
become considerably advective and super critical (Beloborodov 1998). More quantitative 
account of this scenario, however, is beyond the scope of the present paper. 



6. Summary and Conclusion 

We have shown that the multi-color disk blackbody (MCD) emission model can 
successfully reproduce the 0.5-10 keV ASCA spectra of the seven ULXs in nearby spiral 
galaxies. These objects therefore form a homogeneous sample, and are considered to be 
mass-accreting BHs of which the X-ray emission originates from optically-thick accretion 
disks. Their bolometric luminosities, reaching 10 39-40 ergs s -1 even allowing for possible 
distance errors, require the BHs to be relatively massive, up to several tens solar masses. 

The innermost disk temperatures of the ULXs have been found in the range 
T in = 1.1 — 1.8 keV. These values significantly exceed those typically found in Galactic 
and Magellanic BHBs, and are inconsistent with the high BH mass indicated by the high 
ULX luminosities. Similarly high disk temperatures are observed from the two Galactic 
superluminal jet sources, GRO J1655-40 and GRS 1915+105. This problem cannot be 
explained away by changing the disk inclination, the source distance, or the value of £/t 2 , 
within tolerances. Although an optically-thick advective disk (a slim disk) may account 
qualitatively for these observations, we do not find evidence of luminosity-dependent 
deviations from the predictions of the standard disk scinario. 

We hence suggest that the BHs in these objects are spinning rapidly (i.e. Kerr BHs), 
so that the accretion disks can get closer to the BHs and get hotter. In short, a ULX 
may be a Kerr BH with a mass of several tens M Q , in which an optically-thick accretion 
disk is radiating at a near-Eddington luminosity. This provides the first concrete working 
hypothesis with which the ULX phenomenon can be investigated. 
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Fig. 1. — The background-subtracted ASCA spectra of four ULXs, shown without removing 
the instrumental responses. Except for panel (c), the spectra obtained with the SIS 
(higher counts in lower energies) and the GIS (higher counts in higher energies) are fitted 
simultaneously with spectral models. The best-fit model parameters are given in Table 1 
and Table 2, and the fit residuals are presented in bottom of each panel, (a) Spectra of 
IC 342 Source 1, fitted with an absorbed power-law. (b) The same spectra as panel (a), 
fitted with an absorbed MCD model, (c) The SIS spectrum of M81 X-6, together with the 
best-fit absorbed MCD model. Contamination from SN 1993J is also indicated, (d) The 
Dwingeloo 1 X-l spectra and the best-fit absorbed MCD model, (e) Spectra of NGC 1313 
Source B obtained in 1993, shown together with the MCD fit. (f) The same as panel (e), 
but for the 1995 data. 
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Fig. 2. — The deadtime-corrected ASCA GIS spectra of the BH binaries with super-luminal 
jets, fitted with the MCD model or GRAD model. The signal counts are so high that the 
error bars associated with each data point are almost invisible. The fit parameters are 
listed in Table 1, and described in § 4. (a) Spectrum of GRO J1655— 40 obtained on 1995 
August 15-16 (Ueda et al. 1998), fitted with the MCD model plus a power-law. The iron-K 
absorption line is modeled with a notch at 6.95 keV. (b) The same as panel a, but with 
the MCD model replaced by the GRAD model, (c) Spectrum of GRS 1915+105 observed 
on 1995 April 20, fitted with the MCD model. The continuum is modified by photoelectric 
absorption and smeared iron edge. In addition, four atomic absorption lines are included 
(see § 4.2). 
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Fig. 3.— Relation between the bolometric luminosity (Lboi) and the highest color 
temperature (T in ) of optically thick accretion disks around black holes. Various symbols 
indicate measurements with ASCA and other missions. Dash-dotted lines show constant- 
mass grids, calculated via equation (9), assuming the standard accretion disk. Dashed lines 
are those for normalized accretion rates, calculated similarly via equation (11). In addition to 
the data described in the text, plotted are two intensity-sorted data points for IC 342 Source 
1 (Mizuno 2000), three Ginga data points for LMC X-l (Ebisawa 1991), four Ginga data 
points for LMC X-3 (Ebisawa 1991), seven RXTE data points for GRO J1655-40 (Mendez et 
al. 1998), one RXTE data points for GRS 1915+105 (Belloni et al. 1997), and seven Ginga 
data points for GS 2000+25 (Ebisawa 1991). Assumed distances refer to Table 1. For the 
ULXs, i = is assumed. 
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Fig. 4. — Relation between the innermost radius -R; n of an accretion disk (see Table 1) 
and the mass of its central black hole. The BH mass refers to the optical determination 
(M opt ) if available, or the Eddington-limit mass M E otherwise (displayed as lower limits). 
Filled circles, filled squares, and open circles represent ULXs, Galactic jet sources, and other 
"ordinary" BH binaries, respectively. The error bars for Galactic and Magellanic objects 
takes into account typical distance uncertainties. The point for NGC 1313 refers to the 1993 
data. The solid and dotted lines show predictions of a = 1 and a = 1/6 in equation (7), 
respectively, both assuming £k 2 = 1.18. 
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Table 1: Summary of the MCD parameters of ULXs and Galactic/Magellanic BH binaries. a -* 



Source 


distance 


4 h ) 


Jin 


^bol 


T)TT »f C) 

ori iviass ' 


-Kin ' 




(PC) 




v Ke v I 


l iu el &/ b J 


\ 1V1 &) 


(km) 


ULXs 














M33 X-8 e) 


0.72 M 


(0) 


1.15 + 0.03 


O.O 


\ 9 K 


/IO 4- Q 


IC 342 Source 1^ 


4.0 M 


(0) 


1.77 + 0.05 


115 


> 77 


113 + 10 


M81 X-6 9) 


3.6 M 


(0) 


1.48 + 0.08 


39 


^ 91 
> zi 


00 it 


Dwingeloo 1 X-l s ' 


3.0 M 


(0) 


1 on+0.33 


6 1 


> 4 1 


25+2 


NGC 1313 Source B ^ 


4.8 M 


(0) 


1.47 + 0.08 


57 


> 38 


noil? 






(0) 


1.07 + 0.07 


23 


> 15 


147 ± 25 


NGC 4565 off-center *) 


10.4 M 


(0) 


1.39 + 0.08 


117 


> 78 


185±ff 


NGC 4565 center 


10.4 M 


(0) 


1 ^q+0.32 
l.dy„ .23 


47 


> 31 


89ti 


Galactic Transients with Jets 












GRO J1655-40 i) 


3.2 k 


70 




1.3 


~ 7 


19.4l°i 


GRS 1915+105 fe ) 


12.5 k 


70 


1.937 ± 0.004 


12 


> 9 


3o. 9 ;°i 


Standard Galactic and Magellanic Black-Hole Binaries 








Cyg X-l '- m ) 


2.5 k 


30 


0.43 ± 0.01 


0.25 


~ 10 


90 + 18 


LMC X-l m ) 
LMC X-3 m ) 


55 k 
55 k 


60 
66 


n o 9 c;+0.007 

n 07+0.02 
u - y ' -0.01 


3.1 
1.9 


4-10 
5.0-7.2 


86±* 
4°±3 


GS 2000+25 n ) 


3.0 k 


69 


1.10 + 0.01 


5.7 


5-8 


65.31^ 



°) : Based on ASCA measurements except GS 2000+25. 

b ) : The system inclination, assumed to be i = 0° if unknown (in parentheses). 

c ) : For ULXs and GRS 1915+105, the minimum BH mass to satisfy the Eddington limit is 
shown. For the other objects, optically estimated BH masses are quoted. 

d ) : In all cases, (re-)calculated using equation (5) with £ = 0.41 and k = 1.7. 

e ) : Taken from Mizuno (2000), which updates Takano et al. (1994). The photon index of the 
power-law hard tail is fixed at 2.2. 

: Details are given in Table 2. Results are consistent with the report by Okada et al. (1998). 
9) : Details are given in Table 2. 

h ) : The first and second rows refer to the 1993 and 1995 ASCA observations, respectively; see 
also Table 2 and Petre et al. (1994). 

*) : Taken from Mizuno et al. (1999). 

: The same GIS data as used by Ueda et al. (1998), fitted with the MCD model and a 
power-law. See text § 4.1, Fig.2a and Fig.2b. 

fc ) : See text § 4.2 and Fig.2c for detail. 

: The same GIS data as used by Dotani et al. (1997). 

^ : Results from the MCD plus power-law fits. See text § 5.4 for detail. 

n) . -prom Ebisawa (1991), measured with Ginga on 1989 May 3 near the outburst peak and 
fitted with the MCD model plus a power-law. See text § 5.5 for detail. 
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Table 2: Fits to the ASCA spectra of four ULXs. 



Model 


Parameter ^ 


7V H c) 


x 2 /" 




/bol ^ 


IC 342 Source 1 (average) 












Power-law 


1.90 ±0.05 


9.3 ±0.6 


266.5/135 


15.8 




Thermal Brems. 


6.1 ±0.5 


7.5 ± 0.4 


183.2/135 


13.1 




MCD f) 


1.77 ±0.05 


4.7 ±0.3 


137.4/135 


10.7 


12.2 


Unsaturated Compton 


1.42 ±0.09/23 ±2 


6.4 ±0.7 


144.9/134 


11.5 




Broken Power-Law 


1.4 ±0.2/2.4 ±0.3 s) 


6.3 ±0.9 


156.3/134 


12.1 




M81 X-6 












Power-law 


1.91 ±0.08 


4.4 ±0.6 


152.8/82 


5.16 




Thermal Brems. 


5.4 ±0.8 


3.5 ±0.4 


113.6/82 


4.42 




MCD f) 


1.48 ±0.08 


2.1 ±0.3 


83.5/82 


3.69 


4.10 


Unsaturated Compton 


1.1 ±0.1/29 ±4 


2.5 ±0.2 


80.4/81 


3.94 




Broken Power-Law 


1.3 ±0.2/2.8+^ ») 


2.6 ±0.4 


85.1/80 


4.24 




Dwingeloo 1 X-l 












power-law 


1.87 ±0.26 


1 9 ^+ 3 - 2 
iz -°-2.8 


64.4/59 


1.51 




Thermal Brems. 
MCD f) 


7 q+5.3 

' -°-2.4 

i.8o^;| 


10.4 + ^ 


63.9/59 
64.0/59 


1.27 
1.03 


1.14 


NGC 1313 Source B, 1993 July 










power-law 


1.99 ±0.09 


4.0 ±0.7 


174.2/130 


5.43 




Thermal Brems. 


c i +0.9 
°- i -0.7 


2.6 ±0.6 


143.2/130 


4.60 




MCD f) 


1.47 ±0.08 


0.8 ±0.4 


124.5/130 


3.82 


4.18 


NGC 1313 Source B, 1995 November 










power-law 


2.46 ±0.13 


4.4 ±0.8 


110.3/74 


1.51 




Thermal Brems. 


2.9 ±0.4 


2.3 ±0.6 


91.8/74 


1.44 




MCD f) 


1.07 ±0.07 


0.6 ±0.4 


89.8/74 


1.39 


1.66 



a ) : The joint GIS/SIS fits, except the case of M81 X-6 where only the SIS spectrum is used. 

b ) : Photon index for the power-law model, temperature (keV) for the Bremsstrahlung model, 
T; n (keV) for the MCD model, T e (keV)/r cs for the unsaturated Comptonization model, and 
low-energy /high-energy photon indices for the broken power-law. 

c ) : Hydrogen column density for the photoelectric absorption, in units of 10 21 cm~ 2 . 
^ : The 0.5-10 keV flux in 10 -12 ergs s _1 cm -2 , after removing absorption. 

e ) : The bolometric flux in 10 -12 ergs s" 1 cm~ 2 , calculated from the best-fit MCD model. 

: The values of i?j n are given in Table 1. 
a) ■ A break energy is at 3.6 ± 0.5 keV for IC 342 Source 1, and at 3.3 ± 0.5 keV for M81 X-6. 



